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I. INTRODUCTION 



canonical partition function as follows: 



In two previous papers jl], g] a grand canonical parti- 
tion function for non-interacting polymer chains was in- 
troduced in order to compute the chain segment density 
in the context of an ordered monomer ensemble. Here 
we introduce an angular dependence between polymer 
segments in this formalism. In particular, we illustrate 
computations relating to the stiffness of such a chain in 
constraining geometries. Interest in the effects of stiff- 
ness on such polymers extends from biopolymers to liq- 
uid crystalline behavior as well as synthetic stiff and short 
polyamides @ @ g § 0, § |. 

In this paper we present three different models of con- 
fined chains with stiffness. These all involve specific for- 
mulations of the stiffness and nature of the chains in the 
monomer ensemble. In one case analytic expressions for a 
polymer on a cubic lattice can be obtained. For chains in 
a continuum we show that torsional and flexural rigidities 
can be incorporated naturally by the monomer ensemble. 
In a spherical confining region it is possible to solve the 
associated integral equations numerically. 

We use the previously introduced concepts pi |2j. By 
characterizing a bond position by a vector r specify- 
ing its geometrical center and orientation, a bond fu- 
gacity z(r) = (r \z\ r) and interaction weight w(r,r') = 
(r \w\ r'}, i.e. a Boltzmann factor, can be defined. The 
partition function for N bonds 

E N = J [dr 1 dr 2 . ■ - dr N ] (n \z\ n) {r x \w\ r 2 ) . . . 

. . . x {r N ^i \w\ r N ) (r N \z\r N ) . (1.1) 

This is used in order to write the expression for the grand 



S = l + g/l z{wz) N - x 



l + (l\z(I -wz) 1 



(1.2) 



The vector |1) is the vector of ones, resulting in the 
sum/integral over all spatial (and angular) locations. 
One can write for the number density: 



U{r) = E z{r) SzY) 

which gives 

n(r) = — /l (I — zw)^ 1 r\ z(r) (r (I — wz)^ 1 



(1.3) 



By defining tp and ip 

ip(r) = (I — wz)^ 1 1^} , 
tp(r) = (l (I — zw) 1 r^) , 

one can simplify the expression for the density: 

n(r) _ i/j(r)^(r) 
z{r) E 



(1.4) 

(1.5a) 
(1.5b) 

(1.6) 



The grand can onical partition function is calculated from 
equation (1.2) as 
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5=1 + J dr z(r)iP(r). (1.7) 
The average degree of polymerization is given by 

£ = dr n(r). (1.8) 
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The solution of expressions (JL5a) and (1.5b) plays the 
central role in our calculations of the density. 

Whereas in the previous works ^] the physical in- 
terpretation of the vector r represented the location of 
the junction between any two segments of the polymer 
chains, the formalism is identical when a larger degree of 
freedom is represented by a vector of such a kind. 

The different approaches in this paper involve specific 
formulatio ns of the stif fness which are incorporated into 
equations ( 1.5a ) and ( 1.5b ) for the functions ip in or- 
der to determine an expression for the density function 
(PI. Section shows the computation for the density of 
a discrete polymer confined between two parallel plates. 
In the subsequent Section III a general form for the inte- 
gral equations for ip of a chain with bending rigidity and 
torsion is derived. This method is illustrated by numer- 
ical results for the solution of the integral equations for 
a chain in a spherical container. Finally, a simple exam- 
ple for a differential equation formalism is presented in 



such that 



Section IV 



II. POLYMER ON A CUBIC LATTICE 

We make use of a discrete formalism for stiff chains. 
(A similar model was developed in jn^. It was treated 
in a canonical ensemble without any confinement .) In 
addition to associating a position on a cubic lattice, each 
segment also has one of six possible directions along the 
lattice. The bond direction is added to the previous bond 
position to give the next bond position for the chain. 
We map the states for the directions onto real-space unit 
vectors 

(a\ € {(1| ,.-.(6|} <-> {x,y,z, -x, -y,-z} = 

(2.1) 

in a cartoon representation in which all bonds lie along 
coordinate axes and assign to each pair of bonds the 
weight 



w(ri,cri;r 2 ,cr 2 ) = S 



r 2 - 



1 



n - 
if L 



l . 



tn 



1 . 



< a , if t 

b , if t ai • tfj 2 



(7\ ' ^(72 



1 



= -1 



(2.2) 



The first factor in the expression above constrains the 
segments of the polymer; the second factor is responsi- 
ble for the bending energetics. We express all lengths in 
terms of the step-length of the walk. Due to asymmetry 
under exchange of r\ and r 2 of this inte ractio n, two func- 
tions ip, as defined earlier in equations (1.5a) and (1.5b), 
are invoked: 

ip{r,a) = ^1 {I - zwy 1 r,a) (2.3a) 

ip(r,a) = (r,a {I-wz)~ l 1} (2.3b) 



1 = -J2 [ d 3 r (r,a\w T \r',a')z(r',a')ip(r',a') 

a' 

+tp(r,a) (2.4a) 



1 = - 



dV (r, a \w\ r' , a') z(r' , a')ip(r' , a') 



+ip(r,a) 



(2.4b) 



where d 3 r' is here a delta-function measure that converts 
integrals to sums over half-space lattices. 

We investigate this polymer located between two par- 
allel plates located at r z = ±r z = ±N. At the plates the 
boundary conditions require that a segment of the poly- 
mer be oriented in parallel to the plate or perpendicularly 
away from it, but not perpendicularly into it: 



z(r 2 ,cr) = 



0, \r z \>r° z 

0, r z — r° z and a = 3 

0, r z = —r° z and (7 = 6 

zq , otherwise 



(2.5) 



Due to symmetry the functions tp and ip depend only 
on the z-component of position and on a. The x and 
y components are confined to a fixed large length and 
all thermodynamic potentials normaliz ed ap propri ately . 
Furthermore, by comparing equations ( 2.4a| ) and (2.4b) 
for ip and ip with the weight (2.2) substituted, we con- 
clude that 



V ( r z,v) = ip (r z , (<r + 3) mod 6) , 



(2.6) 



with the convention that ?/>(r z ,6) = ip(r z ,0). Under the 
abovementioned conditions we introduce the convenient 
notation: 



tp(r z ,<?) 



ipl ( r z) , 



Symmetry dictates that 

^|| (fz) = 1p\\ (" 

ipt (r z ) = ip { (- 



a = 1,2,4,5 

(7 = 3 
(7 = 6 



■r z ) and 
-r z ). 



(2.7) 



(2.8a) 
(2.8b) 



Since ip is defined in half lattice constants through the 
weight ( |2.2| ) we shall also use the notation ipi(r z ) = ip^, m , 
where m is an integer or half- integer, etc., interchange- 
ably. 

The following two subsections [I A and BlB contain 



the information required to solve the equations (2.4a) 



and ( 2.4b ) for ip and ip. The reader not intere sted in the 
procedure for solution may skip to subsection II C . 



A. System of equations 



By inserting the bending energy and chain position, 
-r° z + 1 < r z < r z — 1, the factor w of equation (2.2) 
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into the condition for ip, equation ( 2.4b| ). the following 
equations are obtained after suitable translations: 

1 = ip\\ (r z ) [1 - z (1 + 2a + b)] - az ipt (r z + §) 

-az ipi (r z - 5) (2.9a) 

1 = ipt (r z - 5) - bz ipi (r z - 5) 

-z tpi (r z + |) - 4az V>|| (r z ) (2.9b) 

1 = (r, + 5) - bzo^ (r z + i) 

-zoV"! ( r z - 5) - 4az ^|| (r*) . (2.9c) 

These equations are valid away from the two plates acting 
as boundaries to the system. Consequently we shall refer 
to calculations relating to the above conditions as those 
pertaining to the "bulk." By using the expressions from 
the above system (equations (2.9a) - (2.9c)) ip\\(r z ) can 
be eliminated, leaving equations expressing ip-\(r z + 1/2) 
&ndipi(r z + l/2) in terms of ipi(r z — 1/2) and ip\{r z — 1/2). 
By defining the column vector 

it is possible to relate functions of tp at different steps by 
ip(r + 1/2) = C -ip(r -1/2) + D (2.11) 
with the matrices 



D 

with 
Ci 



= - 



Cl C*2 

C3 C4 
-«o(l ~b)D x 



l-z (l + 2a + 6) 
z [1 - z (l + 2a + b)+ 4a 2 z ] 
(1 + C 2 ) 
20(1 - 6) 
6 (1 - z (l + 2a + 6)) + 4a 2 z 



C 3 



[l-z (l + 2a + fe) + 4a 2 z ] 

-c 2 

z (l - 6) (1 - C 2 ) 

1 - z (l + 2a + 6) + Aaz Q 
z [1 - 2q(1 + 2a + b) + Aa 2 z 



(2.12a) 
(2.12b) 

(2.12c) 
(2.12d) 

(2.12e) 

(2.12f) 
(2.12g) 

(2.12h) 



As a consequence any "bulk" values of xp can be cal- 
culated given the values of ip at a point (integer and 
half-integer) on the lattice: 



1pn+m 
' l Pn+m+l/2 



+ (C-ty 1 (C" -t)D (2.13a) 

C"t/Wl/2 

+ {C- ty 1 (C n -t)D (2.13b) 



• According to equation (2.9a) ip\\(r) can be com- 
puted with the knowledge of ip(r ± 1/2). 

• The matrix C can be written in terms of C2 as 
follows: 



1+C2 r 
• -C 2 z (l - b)(l - C 2 ) 



(2.14) 



A simple calculation shows that the determinant of 
the matrix is 1, which means that its two eigen- 
values are inverses of one another. 



B. Boundary conditions and solution 

To determine values of xp it is necessary to use the 
equations (2.13a) and ( 2.13b[ ) in conjunction with the 
conditions at the plates confining the polymer. When 
r 2 > N — i it is necessary to refer to the full equations for 



ip ( |2.4rj| ) rather than the bulk values used in the preceding 
subsection. Equations at the upper plate, for example, 
are readily de rived and r eco rded in Appendix |A|. 

Equatio ns (A4b ) and (A5) relate V'TJV-i to V'J.JV-i an d 
equations (A3c) and (A4a) relate ipfN-1/2 1° iplN— 1/2- 



1 = - 



1 = - 



Aaz 



1 - z (l + 2a + b) 
2 ,,2 



(2.15a) 



4a z, 



(i 



bz ipi,N- 



1 -z (l + 2a + 6) 
4a 2 z 2 



l-2o(l + 2a + b) 
Similarly, one can derive for the bottom plate: 

AazQ 



1 



l-z (l + 2a + 6) 



(2.15c) 



Aa 2 z 2 



1 



z (l + 2a- 
Aazo 



+ bz 1p-\~N+l 



(2.15d) 



l-z (l + 2a + 6) 

4a 2 z 2 u \ , 

l-z (l + 2 a + 6) + S *>-™/* 

A solvable system of equations now remains. We know 
from the boundary conditions given above that the two 
component of the column vector xpN-i are not indepen- 
dent and we also know from equation (2.13a) that xpN-i 



is related to xp^N+i- All remaining values for the func- 
tion ip can be determined from the expressions in ap- 
pendix |A| and from the results of the preceding subsec- 
tion. 

The relationship between tpN-i and xp^^+i is 



XpN-1 = C 



2N-2 



Before commencing on further calculations we note: 



ip—N+i + (C — (C 2N ~ 2 — l) D. 

(2.16) 
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Togeth er with the boundary conditions (2.15a) and 
( 2.15b| ) it is straightforward to determine the value of 



i/>at_i from which all other values of ip can be calculated. 



The left and right matrices, L and R, of C diagonalize C: 



LCR 

where LR 
with L 



and A± = — 



A+ 
A_ 

RL = 1 

L\ L 2 
L 3 L4 

1 " 



A 

A- 1 



zb(l-&)(l-C 2 ) + 



1 + C 2 



z (l - b) 



C 2 



4v(* (1 - 6)(1 - Ca)+ 4-*) 



(2.17) 
(2.18) 
(2.19) 



(2.20) 



Equation (2.16) becomes: 



N-l — 



\l N - 2 



A 



2N~2 



A+-1 



A--1 



LD. 



(2.21) 



The solution is: 



L3X1 



A- 2JV+2 L 4 X 1 



X -2N+2 _ l 



x (\- 2N+2 L 3 + \- 2N+2 L 4 X 3 



A- 1 

N -1 



1 



(L 3 £>i + L 4 D 2 



x 2N-2 



L 2 X\ — L\X\ 



(A 

\ 2N - 2 - 1 
A- 1 



2N-2 



L\ -\- \ 2 ^ 2 L 2 X 2 ) 



(L 1 D 1 + L 2 D 2 ) 



\l x X 2 + L 2 - X 4N - 4 {Li + L 2 X 2 ) (L 3 + L A X 2 ) 1 (L 3 X 2 + L 4 )} 



(2.22) 



where we have defined 



x 2 



= 1 



4a zq 



1 - z (l + 2a + b) 



and 



Aa 2 z 2 



1 - Zq(1 + 2a + b) 



bzn 



(2.23) 
(2.24) 



(They are the constant and coefficient in the boundary condition equations.) Standard methods can be used to find 
complete (albeit lengthy) expressions for Li,... ,L±, 



Clearly, the solution ( 2.22| ) is lengthy to write out in full, although it is given explicitly. It is simpler to consider 



limiting expressions. For the purposes of this we shall choose a specific case where b = a and < a < 1, and 
expand to first order in e = a, i.e. the case of an extremely stiff polymer. Functions for this scenario are labeled by a 
superscript "S" . 
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For the stiff case and for n € Z + we have: 

.-n o 



^U-i 



E n 

where E n 



J 

z " 



0(6 2 ) 

1^(^-1) 

c -ty 1 (c n -t)D 
i 



_2N-2 



a?" 1 [ 1 



z - 1 

1 



z 2 n N - 2 -l)+e 



1 



1 - Zq 



4z 
1 - z 

4z 

h e 

1 - z 1-Zo 



zo X—*a 

I—Zq 



27V-2 



0(e 2 



z 2N ~ 2 1 



Z - 1 



1 - *o 



1 , 2z 
i + e i 

1-Zq 1 — Zq 



1-Zo 



. N-l 
+ z 



1 



1-Zq 



l-z 



N-l\ 



(2.25) 

(2.26) 
(2.27) 
(2.28) 

(2.29) 
(2.30) 



Similarly we have the following approximation for the 
floppy ("F") case where a = b = 1: 



1-iZQ 

zo 

1 



-1 



(2.31) 



C. The grand canonical partition function, the 
density, and average degree of polymerization 

The grand canonical partition function, which also fea- 
tures in the expression for the density (1.6), for the dis- 
crete model is given as usual by: 

3 = 1 + ^2^2 z(r,cr)?p(r, a) 

cr { r } 

= l + 2z (4V>||_iv + ^T,-iv) (2-32) 
+z (i, i) • (i - cy 1 { (1 - C 2N ) tf_ w+1/a 
+ (t-c 2N - 1 )^ N+1 



(t-cy 1 (t-c 



2N\ 



- 2N1 



D 



(1 - Cy 1 (l - C™- 1 ) - (2N - ■ £>} 



It is a number dependent on N, zq, a, and b. 

The parameters of the present model can be under- 
stood to give two different types of behavior when con- 
sidering, for example, the parallel and perpendicular ori- 
entations at the center of the two plates. We wish to 
investigate the ratio of the probability that the segments 
in the middle of the plates have a perpendicular orien- 
tation with respect to the plates to the probability that 
they are parallel to the plates using equation (1.6) 



n(0, T) _ V§ 



(2.33) 



The values of -0| (0) and ^(O) at the center are easily 
calculated according to the scheme in Appendix ^ and 
in preceding sections of the paper. 



We calcul ate t he ratio n(0, |)/n(0, ||) for the limiting case where a and b = a 2 are almost equal to using equations 
( f2~29[ ) and ( ^30|) (the "stiff" case labeled by "S"). 



n b (o,T) z o 
n s (0, H) ~ 



JV-l 



l-2 



1-20 1-20 



1 



l-2 



l-2 



1-ZQ l — 2 



1 +^-' 



l-2 



(2.34) 



For chains with a low fugacity and plates which are far apart, this ratio approaches 1, indicating that the chain 
segments are isotropic in the center. For the stiff case it can be seen clearly that for long chains and small plate 
spacing most of the polymer is parallel to the plates. 



The degree of polymerization ( |l.8D can also be evalu- ated 



N 



(2.35) 
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The summation above also contains the non-bulk contri- 
butions at the edges of the system. 



The Potts-type model which has been discussed in the 
present section can also be investigated from the view- 
point of a set of differential equations. The discrete "bulk 
equations" can be converted into differential equations by 
expanding around r z to second order. The coupled set of 
differential equations can be solved by Laplace transfor- 
mation. The need to solve for the roots of a fourth-order 




polynomial for the inverse Laplace transform means that T (a) 

this method does not bring about much of a simplifica- 
tion of the system. 



III. INTEGRAL EQUATION 

To obtain an integral equation for a system with stiff- 
ness one can employ a system of double labeling of succes- 
sive bonds. Thus one can write for the partition function 



L — 'TV 



da da' dbdb' ... (a \w\ a) (a \z\ b) (b \w\ b) 



(3.1) 



where the a,a',b,b', . . . denote successive bonds. The 
scheme is illustrated in Figure |l|. For the bond labeled 
either 3 or 2' in Figure |l| the unit vector for the bend 
can be computed by means of (r 2 — r 3 ) and (r 3 — r^). 
The torsional angle must be computed by taking more 
vectors into account and can be constructed by inves- 
tigating, for example, (ri — 7*2) X (r 2 — r 3 ) in relation 
to (7-3 — ri) x (7-4 — rs). The full torsional and bend- 
ing energy at position 3 in Figure [j] is then expressed 
through a potential V(ri,r 2 ,r 3 ,r i ,r 5 ) = V (f 3 ) which 
can be written as a weight in the chemical potential fac- 
tor z = z(f 3 ) — exp (— j3V (rs)), with r 3 representing 
the supervector centered around position 3. The role of 
w lies in connecting the other points correctly for the pre- 
ceding and succeeding bond-related angles. It does this 
as follows: 

. . . (r\z\ r) (r\w\ r') (r ' \z\ r') . . . 
= ... (r\z\ r) S (ri - r 2 ) S (r 2 - r 3 ) S (r' 3 - r 4 ) 
xS (r' 4 - r 5 ) (f ' \z\ r') ... 

Part (b) of Figure | shows how the preceding 
mathematical prescription of labeling follows when the 
monomers are viewed as having internal structure, with 
appropriate weights for internal conformations. The 
monomers interact so that appropriate parts of the sub- 
structure coincide. This fixes z and w, respectively. We 
remark that there are several possibilities to incorporate 
flexural and bending terms in the integral equation for- 
malism; we shall illustrate one such way. 

The equation for tp (r) becomes 



tp(r) 



(I-wz) 1 1^ 




FIG. 1: The figure depicts the double labeling scheme in 
part (a). This scheme can also bee seen to emerge from 5- 
function Boltzmann factors connecting "monomers" with five- 
segment structures, as depicted in part (b) of the figure. The 
dashed double arrows indicate which parts of the "monomer" 
12345 must be the same physical locations on the "monomer" 
l'2'3'4'5'. In this way the dashed double arrows are delta 
functions of the respective positions. We show one such delta 
function joining points 5 and 4'. 



leading to 

1 = il)(r 1 r 2 r 3 r i r 5 ) 

- J d 3 r' 5 z (r 2 r 3 r 4 r 5 r' 5 ) tp (r 2 r 3 r±r 5 r' 5 ) . (3.3) 

By assuming simple bending without torsional effects, 
the description can be simplified by making use of three 
consecutive position coordinates: 



ip (rir 2 r 3 ) 
and 
$ (rir 2 r 3 ) 



1 



1 



d r' z (r 2 r 3 r') tp (r 2 r 3 r')(3.4a) 



d r' z (r'rir 2 ) tp (r'rir 2 )(3.4b) 



A model which lends itself readily to an iterative nu- 
merical solution is that of the chain of segments of fixed 
length which is confined to a spherical region. A bending 
probability for two adjacent segments labeled 12 and 23, 
with unit vectors h\ 2 and n 23 , can be assigned 



(3.2) 



P ("12, ^23) = P (1 + "12 • n 23 ) . 



(3.5) 
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This causes the forward direction to be favored, with 
p > 0. The function of z = z (t*i, r 2 , r 3 ) in this model 
is to restrict n — r% = ni 2 and r 2 — r 3 = n 23 to unit 
vectors, and to keep the vectors for the spatial locations 
of bonds ri,r 2 , and r 3 from going out of the confines of 
the sphere. Consequently, one can write 

z(r 1 ,r 2 ,r 2 ) = p{l + n 12 - n 23 ) S (\n 12 \ -1)5 (|n 23 | - 1) 
XI? (R - N) 0(R- |r 2 |) i? {R - |r 3 0.6) 

where R is the radius of the confining sphere and i? is the 
Heavyside step function. 

The spherical symmetry and equations ( |3. 6| ) and (3.4a) 
require the following dependence 



where n 2 3,. 

™23- 



i> (rir 2 r 3 ) = -0 (\r 2 \ , n 23 , r ) (3.7) 
is the radial component of the unit vector 



density function at p=0.4 



2.5e-08 



2e-08 



1.5e-08 



le-08 



5e-09 - 



-1 -0.5 0.5 1 
z-component of unit vector 



nz=+l, r=2 
nz=+l, r=25 
nz=+l, r=29 
nz=- 1 , r=2 
nz=- 1 , r=25 
nz=-l, r=29 



FIG. 2: Plot of the density at different radii for R = 30 and 
p = 0.4. Right and left-sloping lines represent forward and 
backward directions, respectively. 



Results 



In a numerical scheme to solv e the equations, it is 
possible to iterate equation ( 3.4a ) at different values of 
the parameters. B y re writing the basic integral equation 
with P, given by (|3.5|), one has 



tp(r,n r ) = 1 + 



dn p(l + n ■ h )zip(r + n, h ). (3. 



In this equation z ensures that the positions of the bonds 
remain within the spherical region. We find (App. |b|) 
that ip can be split into a sum of two parts in successive 
spherical shells, one of which is only a function of the ra- 
dial distance r and another which is directly proportional 
to the radial component of the unit vector n r multiplied 
by a function of r only. Therefore, to integrate from one 
spherical shell to the next we write: 



ip(r,n r ) = <f>i(r) + n r 4> 2 (r). 



(3.9) 



With this manner of splitting the function ip it is possible 
to divide the system into a number of spherical shells for 
each of which a <j>\ and </> 2 have been defined. Equation 
(3.8) for if> can be iterated until the values converge. The 



integration scheme for the different shells is discussed in 
Appendix |b[ 

With the knowledge of ip the value of S can be com- 
puted, and the density expression is: 



directional component is chosen as lying either radially 
outward n z = +1 or radially inward n z = — 1. For cen- 
tral regions of the sphere we see that the straight con- 
figuration is favored and that the angular distribution is 
more-or-less isotropic, i.e. that both unit vectors lying 
inwards-pointing or outwards-pointing is almost equally 
probable. This changes appreciably at the sides, where 
the radially outward density is considerably lower than 
the inwards-facing case. 

In Figure || the density for a both unit vectors tangen- 
tial to the radial direction and the bond being straight 
is plotted as a function of the radius and the probability 
p. The density decreases towards the boundary of the 
system, located at R = 30, and increases with p. Note 
that the chemical potential is built in through p. 

From these graphs a clear picture emerges of a chain 
which is homogeneous in the center of the confining 
sphere and which becomes depleted at the boundaries. 
At these boundaries a tangential orientation of the seg- 
ments is considerably favored above the perpendicular 
(radial) case. 

Figures [| and || indicate aspects of the normalization 
and dependences on orientation. In Figure ^ the depen- 
dence of the degree of polymerization on the combined 
probability and chemical potential p is shown. 



n{fi\2, |r 2 | ,n 23 ) 



—p(l + fl 12 ■ h 2 3)tp(r2, «12,r) 

X^(r 2 ,n 2 3,r). (3.10) 



This five-dimensional quantity can be plotted in a vari- 
ety of manners. In Figure ^| we plot the density at three 
different radii in dependence on the r-component of n 23 
and both unit vectors lying in the same plane. The other 



B. Possible alternative methods for solution 



Another method to solve the integral equation is mak- 
ing use of an expansion in terms of eigenfunctions, and 
using successive substitutions to determine coefficients. 
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Density as function of radius and probability 



density, nl=x, n2=x 




0.1 

Probability 0.3 



10 Radius 



0.4 



FIG. 3: Plot of the density at different radii for R = 30 and 
and various values of p. 



Xi or degree of polymerization as function of p, R=30 

4e+08 | , , , , , , , 1 

3.5e+08 - 

3e+08 - 

2.5e+08 - 

X 2e+08 - 

1.5e+08 - 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 
Probability 

FIG. 4: Plot of the degree of polymerization at a function of 
P- 



Density as function of unit vector z-projections 



density 



2.5e-08 
2e-08 - 

1.5e-08 
le-08 
5e-09 


-5e-09 [ 




FIG. 5: Density for different z-direction projections at p 
0.4. 



For the case of equation ( |3.4a ) a weight 

z{rir 2 r z ) = Aftf [(r 2 - n) • (r 3 - r 2 )] 2 

x (r 2 - rif exp ^--^ (r 2 - r{f 

x (r 3 -r 2 ) 2 exp ^-i (r 3 - r 2 ) 2 ^3.11) 

coul d be introduced. The first factor after d [which, as in 
( |3.6[ ), confines points to an appropriate region] represents 
the bending interaction term of the type cos 2 9 between 
two bonds, while the peaked functions x 2 exp(— x 2 ) set a 
length scale to the segments. J\f is a normalization. 

The weight z and ip can be expressed as a sum of (gen- 
eralized) Hermite polynomials |Ll|] 

z(rir 2 r ?> ) = (H n \r 1 r 2 r ?> ) (3.12) 

■>P(rir2r 3 ) = y m (rir 2 r 3 \H m ) . (3.13) 

m 

These functions c ould be be inserted into the bending 
integral equation ( 3.4a ) and terms compared. 



IV. DIFFERENTIAL EQUATION FORM 

In jl| the density distribution for a polymer confined 
to a spherical cavity was computed by expressing ip in 
equation (1.5a) as 



zip 



w 



- x (^-l) 



(4.1) 



and by making the appropriate choices for the fugacity 
inside and outside the cavity. The inverse Boltzmann 
factor was taken to be the differential operator of 
which w is the Green function. A suitable choice was 
that of a Yukawa form leading to a Helmholtz operator: 



5(r 



_ x Ae- K \ r - r '\ 



(V 2 ~K 2 )w. 



(4.2) 



A solution can be found readily for this system. 

In order to introduce stiffness a segment of the polymer 
chain is now described by the positions of its ends x, x' 
and the orientation of those ends h, n' . Clearly, for a rod 
(x — x'), n and n' are related. The vector \ r) then has 
the dependence \ r) = \ x, #, (f>). By noting that 



1 



da 2 



r 



= 8 {a -a') (4.3) 



a simple multiplicative, bending "Boltzmann" factor is 
introduced: 



w(x, &,4>; x', ■&',<!>') 



Aab 



4 b - x' 



,-a|-#— &'\-b\ct>-4''\-K\x-x'\ 



(4.4) 
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The substitution 



(4.5) 



will separate the variables. For the spherical container 
we need to solve 

(4.6) 



for the outside and 
d 2 



dd 2 



-47raVz -£C 

x 



&2 ) ( V " - k2 ) % 



(4.7) 



for the inside. For the outside the usual radial solution 
is obtained, whereas for the inside of the sphere we have: 



where 



dd 2 

d<j) 2 
cf_ 
dx 2 



-b' 



-4irAa 2 b 



C n 7] 
2,2 



(4.8) 
(4.9) 
(4.10) 

(4.11) 



Here the radial distance and orientation of parts of the 
chain are completely decoupled. This is physically ac- 
ceptable. We impose the fact that the boundary condi- 
tions are cyclic in that £($ + 2tt) — and £((/> + 7r) = 
C(</>). For the polymer confined to the spherical cavity 
the external solution for r only is required. The results 
for r are identical to the results in [|l] : 



r(ip(r) — 1) = 



K 2 - K' 2 

W 2 



(1 + KR) sinh K'r 
K sinh K'R + K 1 cosh K'R 



(4.12) 
,r<R. 



Since the solutions to equations (4.8) and ( J4.9| ) have to 
be periodic or constant, the condition on the constants 
C£ and c<; are that 



fa 2 < and 
+ b 2 < 0. 



(4.13) 
(4.14) 



a variety of problems with more than positional degrees 
of freedom. 

In subsequent work we shall develop this formalism for 
a path integral formulation to include investigation of the 
effects of lateral interactions. 
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APPENDIX A: LATTICE BOUNDARY 
CONDITION AND CENTRAL VALUE 
EQUATIONS 



Here we employ equation (2.4b) for different values of 
r z and a near the boundaries. 

• For r z = r° z + 1, and for a =t,|, and ||, respec- 
tively: 

1 = ^ T (r z ° + 1) (Ala) 

1 = -zoV^ + ^ + l) (Alb) 

1 = V||(r° + 1) (Ale) 



• For r z 

tively: 



,0,1 



, and for a =|,|, and ||, respec- 



,0 , lx 

z~ 2' 



1 = -4az ip\\(r°) -z i>i(r° - §) 



,o , lv 
z ' 2 ' 



(A2a) 

(A2b) 
(A2c) 



1 = -aar ^;(r°)+^||(r° + i) 

For r z — t™, and for a =f, j, and ||, respectively: 

1 = -6z ^(r°)+^T(^ ) (A3a) 
1 = -4a^||(r°-|)-^(r°-l) 

+Air° z ) (A3b) 
1 = -az o -0i(r° - |) 

+ [l-2 (l + 2a + 6)]V||(r°) (A3c) 



With (1.11) Cr) must be negative. 



• For r z 



o _ 1 

z 2 



, and for a =\ and ||, respectively: 



V. CONCLUSIONS 

We have demonstrated how the formalism accommo- 
dates potentials with angular dependence in three differ- 
ent ways. The formalism is quite generally applicable to 



-Aaz ip\\(r° z ) 
-bz ipi(r" z - 



1 = [l-zb(l + 2a + &)]tf||(r°-£) 
-azoi>i(r° - 1) 



(A4a) 



(A4b) 
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R-2 R-l R 




FIG. 6: Schematic representation of the integration proce- 
dure over shells. The shell numbering is illustrated, as well 
as the mixing between different shells. "F" represents the 
forbidden region for any bond of the polymer. 



For r z = r® — 1 and a =\: 

1 = -4az ip\\(r° z - |) 

-bzMr° z -l)+^{r° z -l) (A5) 



Halfway between the plates the symmetry dictates that 
V^o = ipi,o- (A6) 



APPENDIX B: NUMERICAL INTEGRATION 
SCHEME FOR SPHERE-CONSTRAINED WALK 



In order to elucidate the integration scheme used for 
the numerical calculations we refer to Figure | The 
spherical geometry of the system is shown here up to 
the edge of the system. For each shell of thickness 1 we 
calculate ipi = 4>f^ + n z <pp . Since the shells have the 
thickness of the radius of the bond, at each stage there 
are contribution s fr om the two adjoining shells, accord- 
s' These can be summed if we assume 
f>^ and <jj?' are constant and approxi- 



ing to equation ( 3 



that the values of < 
matcly equal to the value of the </>'s in the middle of each 
shell. In the final shell ending at the radius R, the bond 
vector is permitted to move only in the allowed region A 
with zero weight in the forbidden region F of Figure ||. 

The set of , <f>f^ j is iterated through equation (|3.8| ) 

until the values no longer change. 
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